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'. ABSTRACT 

Context. Very few examples of luminous blue variable (LBV) stars or LBV candidates (LB Vc) are known, particularly at metallicities 
below the SMC. The LBV phase is crucial for the evolution of massive stars, and its behavior with metallicity is poorly known. V39 in 
fT^ \ IC 1613 is a well-known photometric variable, with B-band changes larger than Imag. over its period. The star, previously proposed 

to be a projection of a Galactic W Virginis and an IC 1613 red supergiant, shows features that render it a possible LBVc. 
Aims. We aim to explore the nature of V39 and estimate its physical parameters. 

Methods. We investigate mid-resolution blue and red VLT-VIMOS spectra of V39, covering a time span of 40 days, and perform a 
] quantitative analysis of the combined spectrum using the model atmosphere code CMFGEN. 

• . Results. We identify strong Balmer and Fen P-Cygni profiles, and a hybrid spectrum resembling a B-A supergiant in the blue and 

i-^ ' a G-star in the red. No significant Vrad variations are detected, and the spectral changes are small over the photometric period. Our 

analysis places V39 in the low-luminosity part of the LBV and LBVc region, but it is also consistent with a sgB[e] star. From this 
analysis and the data in the literature we find evidence that the [ff/Fe] ratio in IC 1613 is slightly lower than solar. 
, Conclusions. The radial velocity indicates that V39 belongs to IC 1613. The lack of Vrad changes and spectroscopic variations 

^ ■ excludes binary scenarios. The features observed are not consistent with a W Virginis star, and this possibility is also discarded. We 

^ ' propose that the star is a B-A LBVc or sgB[e] star surrounded by a thick disk precessing around it. If confirmed, V39 would be 

the lowest metallicity resolved LBV candidate known to date. Alternatively, it could represent a new transient phase of massive star 
evolution, an LBV impostor. 
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1. Introduction sity, going from intense and sudden eruptions in which the stellar 

magnitude increases by two magnitudes or more (like rj Car) to 
^ i Luminous blue variables tl.RVs) constitute a short and rare i- • i ^ •. ^: ^_ 



X 



Lummous blue variables (LBVs) constitute a short and rare ^y^jj^ variations of 1-2 magnitudes in time scales of years to 
^ phase in the life of massive stars, in which the stars suffer ^^^^^^^ (jjj^g ^he prototype S-Dor). These photometric changes 
^ great mass loss with unpredictable outbursts and mass eruptions. accompanied by the corresponding spectroscopic ones, in 
Apart from the stages in which a nuclear runaway is produced in ^^1^^ effective temperature varies and the bolometric lumi- 
the stellar core, LBVs seem to be the shortest phase in the mas- j^^sity may or may not remain constant ( Vink 2009). However, 
sive star's evolution, even shorter than Wolf-Rayet phases, with outbursts and eruptions are not very common, and sometimes the 
_ estimated hfetimes of 10^-10^ years. In spite of its brevity, this .^^^ ^^^^^^^ ^^^^le for centuries, as it is the case for P-Cygni. 
d • phase IS extremely important for the evolution of the star and ^hus, rather than from photometi'ic valuations, they are often 
■ ■ ' Its surroundings, because it returns important amounts of pro- searched for by means of their P-Cygni profiles in the Balmer 
cessed material to the ISM, often forming a shell or a nebula jj^g^ ^^^ich are a signature of the presence of strong mass- 
around the star. Moreover, the loss of matter modifies the star's j^^^ processes), emission line profiles (at lower resolutions or 
subsequent evolution and it has been very recently suggested ^IR wavelengths) and of other lines which depend on the tem- 
thatLBVs might be direct SNe progenitors (see | Kotak&Vink | pgrature of the star When the star shows a chai-acteristic spec- 
(l2006!),[VQiikJ (12009D), which may have an impact on our current ^.^^^ ^j^^ p.Cygni profiles, but no relatively large photomet- 
ideas of massive star and galaxy evolution in the Universe. spectroscopic variations (or we do not have information 
It IS a tautology to state that LBVs are luminous (log(L/Lo > about such variations), we speak of an LBV candidate, or LBVc. 
5.4), blue ((B-V)o < 0.1) and variable stars. Nevertheless, their Sometimes, other signatures (like IR excess or nebulosity, see 
variations can be considerably different in frequency and inten- f.e. Clark et al. ( 200^)) are considered further signs for the clas- 
sification as LBVc. 
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* Based on observations obtained at the ESO VLT for Programmes The mechanism producing LBVs is not known, neither is the 

078.D-0767 and 080.D-0423. exact duration of the phase (and hence, the exact impact on the 
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Star's life). It is assumed that the stellar mass and, very espe- 
cially, the mass-loss rate play a role. As the winds of massive 
stars are driven by the momentum gained by the atoms from the 
radiation field through the absorption of photons by metal lines, 
LBVs should be rare at low metallicities. This is consistent with 
no LBVs or LBV candidates known in the Local Group at metal- 
licities below that of the Small Magellanic Cloud (SMC). While 
this fact could be at least partly the result of the low number of 
massive stars in such galaxies, the sheer presence of an LBV at 
low Z should be considered significant, as their presence offers 
clues for mass-loss processes that could have been relevant in 
the early Universe. 

Beyond the Local Group few r elevant examp l es of 
LBVs or possible LBVs ar e known. iDrissen et al. J (1200 Ih 
(see also iPetit et al~l (12006')) have analyzed the star VI 
in NG C 2363, a galaxy a t abou t 3.3 Mpc (iTolstov et al. I 
1995h. [K arachen tsev et al. I (l2002h ) with 12- Hlog(0/H)^ 7.9 



Gonzalez -Delgadoet al. II 19941) or Z= 0.25 Zq dLuridiana et al. I 
999h . HST spectra were used for the analysis to avoid con- 



tamination because the star is in a dense stellar environment 
with strong nebulosity. NGC 2 363- VI shows str ong photomet- 
ric and spectroscopic variations (iPetit et al. l2006l) . thus fulfilling 
all requirements to be called an LBV. Us ing the CMFGEN code 
(Hillier& Miller 1998), Drisse n et alTI (lOOl) obtained Teff= 
11000 - 13000 K, log(L/L0)= 6.35 - 6-40 and Voo= 325 - 290 
km s ' for spectra collected in Nov. 1997 and July 1999, re- 
spectively, fully in agreement with the values expected for an 
LBV in transition. From the iron spectrum they get 0.2 
Zpi, intermediate betw een Gonzalez-Delgado et al._ ( 1994) and 
iLuridiana etsHim^ . 

Very interesting, but less evident is the case of two very 
recent objects studied in very metal po or galaxies: knot 3 
in DDO 68 (Pustilnik et al. (2008); Izo tov & Thuan I (|2009|) ) 
and PHL 293B (Izotov & Thuan . ,2009.) . DDO 68 is a blue 
comp act dwarf (BCD) galaxy at a dis tance of at least 5.9 
Mpc dMakarova & Kar achentsev I (Il998l) determined its dis- 
tance from the average magnitude of the three brightest blue 
stars). K not 3 (DDO 68 - 3) is one of its strong HII regions, 
for which'Pustilnik et al. ' OOO^ give 12H-log(0/H)= 7.10+0.06 
and Izotov & Thuan (2009) 7. 15+0.04. In spe c tra of this knot 
separated by about three years, Pustilnik et al. (2008) detected 
spectral variations (broader emission in the Balmer lines and 
He I emission not previously present) at resolutions of 6 and 12 
A FWHM that they attributed to the outburst of an LBV in- 
side the region. This explanation was adopted and further inves- 
tigated by llzotov & Thuan I (l2009t) . who estimated the absolute 
magnitude, H„ luminosity and wind terminal velocity from APO 
and MMT spectra at resolutions of 7 and 3 A FWHM. How ever, 
note that the 1.5" slit width used by Izotov & Thuan (20()9) cor- 
responds to about 45 pc at the distance of DDO 68 and therefore 
the observed spectrum is a composite of the HII region and the 
underlying stellar population (including the possible LBV). 

A simi lar problem is f ound in the second object pre- 
sented by llzotov & Thuan! ( |2009|) . PHL 293B. This is an- 
other low-metalhcity B CD ( 12+log(0/H)^ 7.66+0.04 ac- 
cording to [izoto v et ak] (|2007|) and 7.72+0.01 according to 
llzotov & Thum I (12009)) bu t this time locat ed at 2 1 .4 Mpc (NED 
database). The spectra of Izotov & Thuanl (2009) were obtained 
with a slit width of 1.0", which corres ponds to more than 100 
pc. The oxygen abundance quoted by llzotov & Thuan I (|2009|) 
has been obtained from a high resolution (0.2 A FWHM) VLT- 
UVES spectrum that also displays broad emissions at the po- 
sitions of H„, H^ and Hy and blue-shifted absorptions in all 



Balmer lines. No He broad emissions are seen. Neither the DDO 
68-3 nor the PHL 293B spectra show metal lines that could be 
attributed to a star. In both cases, the wind terminal velocities 
quoted bv llzotov & Thuan I (l2009l) (» 700-850 km s ') are sig- 
nificantly larger than the usual range of LBVs ( 100-250 km 
500 km s"' for 77 Car), which is attributed by the authors to the 
very low metallicity. 

The mere presence of LBVs at these low metallicities rep- 
resents a challenge for the theory, as there are few metals to 
drive the wind. Other processes have been suggested to cause 
the LBV outburst, like the metallicity-independ ent continuum- 
driven (instead o f line-driven) wind proposed by lOwocki et al. I 
(2004) (see also ISmith & Owocki I a006) ). Note however that 
the star needs to be close to the Eddington limit for the mecha- 
nism to operate. There is a clear interest in investigating these 
processes, confront them with observations and apply the re- 
sults to early Universe and Population III objects. Obviously, it 
would be very interesting to obtain spatially resolved spectra of 
the stars in DDO 68-3 and PHL 293B, both to confirm that the 
changes in the composite spectra are due to an LBV and to study 
the evolutionary processes at very low metallicity. However, be- 
cause of the large distances involved, we have to look closer for 
resolved stars, such as NGC 2363-Vl. 

IC 1613 is a dwarf irregular galaxy in the Local Group wit h 
a distance modulus of (m-M)o= 24.27 (Dolp hin et"ani2001h . 
It has a metallicity of log(0/H)+ 12= 7 80+0 .10 as determined 
from its B-supergiants (Bresolin et al. I 2007b. whereas nebular 



studies vary between 7.60 and 7.90 (iLee et 311120031) . It is there- 
fore clearly below the SMC metallicity and intermediate be- 
tween that of NGC 2363 and PHL 293B. Therefore, the obser- 
vation of a resolved LBV in such a poor-metal galaxy would 
constitute a step forward, adding a second object to the list of 
LBVs at metallicities below that of the SMC and lowering the 
present metallicity limit for resolved LBVs. 

IC 1613 shows a recent and intense burst of massive star for- 
mation, particularly in its NE part. We have recently published a 
new catalog of OB associations in IC 1613 dGarcia et al. II2009I) 
and their physical properties (Garcia et al., in prep.) as part of 
our effort to carry out an in-depth study of the young population 
of IC 1613. We also obtained spectra of some stars in this galaxy. 
In the field of IC 1613 we find the well-known variable star V39. 
In Table [T] we offer an overview of its photometric data. 

The nature of V39 has remained unclear since its discovery. 
Sandage ( 1971), in his analysis of IC1613 based on unpublished 
previous work by Baade, points out that V39 is the only peculiar 
variable in the field, because of its inverted yS-Lyrae light curve. 
This was already noticed by Baade, who never considered it to 
be a Cepheid, not only because of this peculiar light curve, but 
also because it was too bright to fit the P-L relation compared to 
other Cepheids with the s a me pe riod. Consequently, in spite of 
a suggestion by ISandagel (Il97ll) (that implied a new P-L rela- 
tion and a change in the distance modulus to bring it into agree- 
ment with other Cepheids) mos t auth ors did not consider it a 
genuine Cepheid. lUdalski et al. I (1200 ll) used its position in the I 
vs (V-I) color-magnitude dia gram as a new argument to discard 
it as a Cepheid. Very recentlv JPietrzvhski et al. I(l2006 l) have also 
shown that V39 is too bright to fit the Cepheids P-L relation in 

the infrared. 

lAntoneUo et al. I d 19991) showed that the combined light 
curve from different observers is the result of the superpositon 
of two periods: a long period of 1123 days and a short period 
of 28.699 days. They showed that by subtracting the long pe- 
riod from the light curve, the remaining one did not show a 
clear difference between the primary and secondary maxima, so 
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Table 1. Photometric data for V39. IUdalski et al. I (l200l . |Pietrzvriski et al. I (l2006h and lGarcia et al. I (l2009h give mean values. 
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that the short period cou ld be halved. Based on these results, 
iMantegazza et al. I(l2002h proposed V39 to be actually the casual 
overlapping of a distant Galactic W Vir star and a red supergiant 
belonging to IC 1613. If confirmed, the W Vir star would be at a 
distance of at least 1 15 Kpc, being the most distant Galactic star 
from the Galactic plane, and raising the question about the true 
extent of the Galactic halo. 

In this article we show that V39 does actually belong to 
IC 1613 and analyze its possible nature, including the possibility 
that it is an LBV candidate, the first one known in IC 1613. If 
confirmed, it would be the (resolved) LBV candidate within the 
lowest metallicity environment. The observations are presented 
in Sect.|2] a description of the spectrum is offered in Sect.|3]and 
the nature of V39 is then discussed in Sect|4] An estimation of 
the stellar parameters is offered in Sect. |5] and the evolutionary 
status and other properties are discussed in Sect. |6l Finally, the 
conclusions are presented in Sect. |7] 

2. Observations and data reduction 

Observations were performed with VIMOS at VLT in MOS 
(multi-object spectroscopy) mode. The HR-Blue and HR- 
Orange gratings were used, resulting in a resolution of R~ 2050 
and R~ 2150, respectively, and a spectral wavelength coverage 
of 3870 to 7240 A. In total 3x19 blue spectra were observed 
from October 5 to November 5, 2007, and 2 x 10 red spectra 
from November 5 to November 13, 2007. Integration times for 
individual exposures were 650 sec for the blue spectra and 1005 
to 1 1 10 sec for the red spectra. The three (two) consecutive expo- 
sures in the blue (red) were coadded. The coadded spectra were 
extracted with standard IRAlQ procedures after wavelength cal- 
ibration and cosmic ray removal. Next, a barycentric correction 
and a c orrection for the systemic velocity of IC1613 (-234+1 
km s"'. lLuetaI.1(ll993h ) were applied. Finally, all 19 (10) re- 
sulting spectra were coadded to a single blue (red) spectrum. 
The resulting spectra, with an average SNR of 80, were then 
rectified. 

Figure [T] shows the finding chart of V39. The star is located 
at Q'(2000)= Olh 05m 2.04s, 5(2000)= 02'''=g 10' 24.7" (from the 
astrometry by Garcia et al. ( 2009)). It clearly lies in a zone of in- 
tense recent star formation, with strong bubbles visible in Fig.[T] 
Table [T] gives the main photometric data of the star, collected 
from the literature. Details for our astrometry and photometry 
are given in Garcia et al. (2009). At the time of the photometric 
observations (July, 30th, 2006), the star was in a state of min- 
imum brightness (see Sect. |4|. Note that photometry and spec- 
troscopy were obtained with a time difference of more than one 
year. 

Unfortunately, no HST, 2MASS or Spitzer images of V39 
are available, and Galex does not detect V39 in either the FUV 

' IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 



or the NUV bands. Thus the data available for this work extend 
from the U to the K band. 



3. Spectral morphology 

Figures |2] and |3] show the observed coadded blue and red spectra 
of V39, respectively. The absence of [O iii] or He i emission lines 
indicates that there are no contamination problems with nebular 
lines. The spectrum is dominated by the strong P-Cygni profiles 
of the B aimer lines. They are shown in detail in Fig. [H where 
we can see the strong increase in the emission profile through 
the series. The second outstanding characteristic is the relatively 
strong P-Cygni profiles of the Feii lines at AA 4924, 5018 and 
5168. The first two are very close to the Hei lines at AA 4921, 
5016. However, the fact that other He i lines are much weaker or 
even absent, together with the presence of Fe ii /15 168, indicates 
that these are Fe ii lines. The only He i line seen in the spectrum is 
He I /15876. This line is present in the red and blue spectrograms, 
which overlap in a wide wavelength region. In Fig. |5] we plot 
part of this overlapping region, which shows a weak P-Cygni 
profile of He i /15876 in both spectra. The presence of this He i 
line is very important as it can be used to constrain the stellar 
temperature and He content. 

As can also be seen in Fig.|5] many weak features are present 
in both spectra. This shows that these features are real and not 
the product of a noisy spectrum. The strong absorption lines red- 
wards from He i are the Nai D lines at AA5890, 5896. The fact 
that they fall exactly at their laboratory position after barycentric 
and systemic velocity corrections indicates that they belong to IC 
1613, not to the Galaxy and that they are at rest with respect to 
the star. We also notice a weak redshifted emission component 
for both lines, which can initially be identified either as the emis- 
sion component of a P-Cygni profile or as the corresponding Na i 
sky line at the observer's rest frequency (therefore displaced by 
the IC1613 systemic velocity in the stellar spectrum). Checking 
individual exposures, we realized that in two of them we could 
not completely remove the sky lines, and they appear as weak 
redshifted emission components in the stellar spectrum. 

The third characteristic is that the observed spectral type 
seems to correspond to later types when one moves from blue to 
red spectral ranges. This is shown in Fig.|6l where we compare 
V39 with an A 3I star and a G9I star , extracted from the UVES 
POP database ( Bagnulo et al. 2003) and degraded to VIMOS 
resolution. Apart from the difference in line intensities (which 
can be partly attributed to the difference in metallicity) we see 
that V39 follows more the A3 spectrum in the blue range and 
the G9 in the red one. Note that this is not an attempt to classify 
the star, but just an illustration of the spectral differences. 

As V39 is a photometric variable, we would expect the spec- 
trum to also vary periodically. However, changes in our individ- 
ual spectra are small and mostly compatible with difficulties in 
the sky line reduction or within noise limits. For example, the 
variability we can see in the Fe ii profiles is modest. Moreover, 
no significant radial velocity variability is detected (although 
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Fig. 1. Finding chart for V39: R band VLT-VIMOS image of IC 1613, showing the active star-forming NE lobe of the galaxy 
(left) and a zoom into V39's neighborhood (right), indicated by the square in the left panel. In both panels the position of V39 is 
marked with a cross. The slit for the observation was placed in the North-South direction; there is no resolved nearby star that could 
contaminate the spectrum in this direction. 




Fig. 4. Balmer lines (H^ to Ha) in the spectrum of V39. Note the 
scale change in the y-axis for Hq,. Wavelength is given in A. 

the relatively low resolution and modest S/N of the individual 
spectra limit our ability to detect radial velocity changes). As 
an example we can examine the line profiles of H„, as this line 
is known to be very sensitive to conditions in the atmospheric 
layers. In Fig. [T] we plot the profiles observed during seven 
days (i.e. covering a quarter of the 28 days period or half the 14 
days period), and we see that only small changes are present. H/?, 
observed during a time span of 30 days, displays similar small 
changes. Figure [8] displays the equivalent widths measured in 
the individual H,,. spectra for reference. We see again that the 
variations are small. Changes in are once more similar, al- 
though because of the smaller values and lower SNR, the relative 
changes appear a bit larger. 

4. The nature of V39 

The spectrum described in the previous section resembles partly 
that of an LBV candidate. The strong emission and moderate 
maximum velocity of the P-Cygni profiles of the Balmer lines 
point to an intense, dense and slow wind, like that of a super- 
giant. At the same time, the Fe ii lines and the lack of a strong 
He spectrum points to an intermediate temperature. The cool fea- 
tures in the red spectrum indicate some degree of contamination 
by a secondary source that however shares the radial velocity of 
the blue spectrum. In addition, V39 displays periodic or quasi- 
periodic photometric variations that are not expected in an LBV. 
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Fig. 5. Comparison of a region of lC1613-LBVc-l observed in 
the red (upper) and blue (lower) spectrograms. The He i and Na i 
line positions have been marked. 



For that reason we analyze here different alternative scenarios 
for V39 using the previously known facts about V39 given in 
the introduction and the spectral description from previous sec- 
tion. 

i) a Galactic W Yir star plus an IC 1613 red supergiant . 
Following Antonello's et al. ideas, Man tegazza et al. I (l2002h 
obtained two periods of 14.341 and 1118 days. They offered 
a solution for V39: the superpositon along the line of sight of 
a Galactic W Vir star, a Population 11 Cepheid pulsating with 
a period of 14.341 days, and a red supergiant in IC 1613 that 
varies with a period of 1118 days. The alignment has to be 
in a way that both stars contribute to the light we receive at 
Earth. This would agree with a composite spectrum and, in 
spite of the unprobable configuration, this interpretation has 
an additional interest: it implies that the W Vir star should 
be at a distance of 115 Kpc. This would mean that it is the 
most distant Galactic star known and would raise the ques- 
tion about the true extension of the Galactic halo. 
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Fig. 2. Blue spectrum of V39. The P-Cygni profiles of Hy, H^, and Fell] 4924, 5018 can be easily recognized. Note the absence of 
Hel lines. Wavelength is given in A. 
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Fig. 3. Red spectrum of V39. Wavelength is given in A. 





Fig. 6. Comparison of the blue (above) and red (below) spectral ranges of V39 (black solid line) with the spectra of HD 62623 (A3I, 
green) and HD 1 17440 (G9I, red, dashed). 
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Fig. 7. Time series of Hq, spectra of V39. Relative fluxes have 
been displaced according to the time difference with respect to 
the first observed spectrum. Units in the ordinate axis are half- 
days. The time difference between the first and last spectra is 
seven days. Wavelength is given in A. 

However, our higher resolution spectrum is not consistent 
with this classification. Firstly, the barycentric plus systemic 
velocity correction is consistent with the rest of the stars ob- 
served in IC1613, both in the blue and red parts of the spec- 
trum, and therefore consistent with V39 being a member of 
IC 1613. Moreover, there is no radial velocity shift between 
the red and blue spectra. And secondly, the timeseries of the 
individual spectra before coadding them shows that the P- 
Cygni profiles of the Balmer lines vary only slightly (see 
Fig. |7]i, whereas Balmer lines in W Vir stars vary strongly 
alon g their cycles and neve r show P Cygni profiles (see Fig. 
3 of iLebre & GiUetl (Il992h ^. 

ii) A close binary system. Another possibility to explain the pe- 
riodic variations is V39 being actually a binary system, com- 
posed for example of a supergiant dominating the spectrum 
and a hotter companion. We could assume the binary system 
motion being responsible for the short period, whereas the 
supergiant would be responsible for the long-period variabil- 
ity. However, this explanation is not supported by the data. 
The short period of about 14 or 28 days is too short to be 
caused by binary motion. It would imply separations smaller 
than 75 or 150 Rq respectively (which are not compatible 
with the presence of a red supergiant) and large radial veloc- 
ity variations that are not detected. 

iii) A wide binary system. Alternatively, we could attribute the 
long period with its small magnitude changes to the bi- 
nary motion, and the short period with its large magnitude 
changes to the primary (brightest) star The long period of 
1120 days would imply a larger separation than the short 
one, about 1800 R©, and modest radial velocity variations 
(of 30-50 km s"', which means a maximum separation of 
about 60-100 km s"' every 1.5 years). At the distance of IC 
1613, this spatial separation is too small to be detected, and 
our spectral resolution is too low to detect the implied ve- 
locity changes in only 40 days (that is the time span covered 
by our observations). Assuming that the light curve of the 



long period is due to binary motion, we would attribute a 
change of AB= 0.3 to contamination by the secondary (from 
Figure 1 1 of Antonello et al. ( 1999)), corresponding to a star 
with mg > 20 or Mg > -4.2, which is too low for any super- 
giant, implying an object with a small radius and a blue color. 
Depending on the exact relative magnitudes, this may or may 
not be detectable in the composite spectrum. However, the 
major photometric changes should come from the supergiant 
and are thus inconsistent with no spectroscopic changes dur- 
ing the photometric period. 

iv) Contamination by a nearby object. This situation is similar 
to that of a binary system, but without the restrictions on ob- 
ject separation. Therefore, we will only consider the case that 
the short period variations are due to a nearby star or binary 
system, close enough along the line of sight to the star that 
dominates the observed spectrum (i.e., much in the line of 
Mantegazza et al., but now with both objects belonging to IC 
1613). In this case the spectrum would be dominated by a red 
supergiant responsible for the long period and the second ob- 
ject would be a star or stellar system in IC 16 13 responsible 
for the short period variations, either because it is an eclips- 
ing binary system or because it is a pulsating star. However, 
as the magnitude changes are considerable, the nearby ob- 
ject should be relatively bright, which restricts the possibili- 
ties. If, for example, we attribute the AB= 1.0 mag. and the 
AV = 0.5 variations to this n earby object (using Figs . 11 an d 
3 of I Antonello et al. I d 19991) and [Mantegazza et al. I ( |2002|) . 
who have already substracted the long-period effect) the ob- 
served photometric changes are compatible with the super- 
giant dominating the blue spectrum only during less than 
half of the period. As we have observed our object for more 
than one period (when we take blue and red spectra into ac- 
count) without significant spectral changes, we discard this 
scenario. 

v) A star surrounded by a hot thick disk. In this scenario the 
star is surrounded by a thick disk that emits at a characteristic 
temperature of 5000 K, producing the spectral features in the 
red part of the observed spectrum. As we do not expect a star 
of about 10000 K to show relevant spectral features in the 
6000 A region at our resolution (see Fig.|6l), we can attribute 
them to the disk even if it does not completely dominate the 
system light at these wavelengths (of course, this is not true 
for some particularly intense lines in the spectrum of the star, 
like Hq,; conversely, the stellar Balmer lines in the blue could 
also be affected by the disk lines). The blue spectrum would 
be produced by the star itself and represent a hotter tempera- 
ture. The photometric changes would be due to disk preces- 
sion (at least partly) changing the aspect angle of the sys- 
tem. The relative contributions of star and disk would vary 
periodically, but always producing the same combined spec- 
trum, except for small spectral line variations (in absence of 
precession, photometric changes without spectroscopic ones 
would imply a change in the stellar bolometric luminosity, 
which also argues against a single star explanation). In this 
scenario, photometric changes are wavelength dependent, as 
observed, as the disk becomes thinner at larger wavelengths. 
Details of such a scenario have still to be worked out to show 
that the large photometric changes are actually compatible 
with the small spectroscopic changes, but it is consistent with 
present data. Moreover, the possible presence of a disk is re- 
inforced by the large reddening of V39, much larger than 
the average one towards IC 1613, than the typical internal 
reddening in this gala xy, and than that of the surroundings 
of V39. For example, iLee et al. I (fT993i) give E(B-V)=0.02 
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2454410.0 (DAYS) 



Fig. 8. The equivalent widths of Ha (in A) of the individual spec- 
tra available. Upper plot: absorption profile; lower plot: emission 
profile. Changes do not show any clear periodicity, nor are they 
larger than a few percent. 



towards IC1613, while iGarcia et al. I (|2009|) give E(B-V)= 
0.10-0.15 as the mode of the distribution of extinction for 
OB associations. Moreover, the spectra of the three stars 
closest to V39 in our VIMOS mask do not show apprecia- 
ble Nai lines. Therefore the reddening must be very local. 
However, we have no further independent evidence support- 
ing the disk scenario. New spectra with higher spectral res- 
olution to investigate for example the presence of a double 
peak produced by a disc would be very helpful to confirm or 
discard it. They should be ideally obtained simultaneously 
with new photometric data, as the available photometry has 
been obtained in the course of decades, by different authors 
using different filters. 



5. Spectroscopic analysis 

We have concluded that the observed spectrum of V39 is best 
interpreted as that of a BA-type supergiant surrounded by a hot, 
thick disk precessing around it. The star itself could be an LBV 
candidate, attending to the strong P-Cygni profiles shown in the 
Balmer and Fe ii lines. To be classified as an LBV the star should 
show outbursts and photometric and spectroscopic variations. 
There are no known outbursts of V39, but part of the periodic 
or quasi periodic photometric variations co uld be intrinsic to th e 
star This possibility has been explored by Onifer et al. ' ('2008'), 
who show that it is possible for LB Vs at low metallicity and large 
He abundances to develop quasi-stable pulsational instabilities 
with periods of about one month due to a A--effect mechanism 
caused near the iron bump opacity. 

To further investigate the possibility that the star is an LB Vc 
we have carried out a quantitative analysis of the spectrum to 
place it in the theoretical Hertzsprung-Russell Diagram and to 
compare its position with those of other stars. For this purpose, 
we performed a model atmosphere analysis of the spectrum of 
V39 obtained after co-adding the individual spectra, as described 
in Sect.|2](19 blue spectra and 10 red spectra). 

The stellar parameters were deter mined by means of th e 
CMFGEN model atmosphere code (iHillier & Miller I Il998h . 
Details of the versi on and atomic models used are given in 
iNaiarro et al. ] ( I2OO9I) . In Tables |2] and [3] we list the main model 
parameters and the synthetic photometry, obtained by integrat- 



ing the emergent spectral distribution weighted by appropriate 
filter values. 

As no ionization equilibrium could be used for Tefj determi- 
nation, we have adopted a He abundance charac t eristic of LBV 
stars, namely H/He= 2.5 ( iNaiarro et all d 19971) . iNaiarro et aP 
(I2009h 1. Teff was then determined paying attention to the pres- 
ence of Fe 11 P-Cygni lines, the probable presence of He i /15876 
and the absence of other Hei lines. We choose to fit the H(j, line 
as this is the strongest spectral feature, and should be dominated 
by the mass-losing star in spite of the contamination of the red 
spectrum. 

The radius (and therefore the luminosity) was derived from 
the K-band magnitude, as this is the wavelength with the lower 
reddening and for which the magnitude varies less (we note that 
the changes in brightness decrease with increasing wavelength). 
We assume here that in a star4-disk scenario the stellar contribu- 
tion dominates at this wavelength. To o btain a reddening es tima- 
tion we used the photometry given by lUdalski et al. I (12001.) . as 
it covers a time span similar to that of our spectroscopic obser- 
vations (40 nights in our case, 51 in that of Udalski et al.; note 
that using our own photometry given in Garcia et al. (2009) is 
less appropriate, as it was obtained in one single night). 

We proceeded in the following way: from the unreddened 
(V-I) color obtained f rom the model and the observed (V-I) from 
lUdalski et all (12001') we derived (Av - Ai)= 0.92. Using then 
(Ai/Av)= 0.601 and (Ak/Av)= 0.112, we finally got Ak= 0.26 
and K= K^odei + 0.26= 15.73 (compared with the observed 
K= 15.82). With these values we obtained 1= 17.30 and V= 
18.53, which i s sli ghtly brighter than the values observed by 
lUdalski et al.1 12001), namely 1= 17.62 and V= 18.85. This is 
more clearly seen in Fig.|9l where we plotted the observed mag- 
nitudes from Table [1] and the predicted magnitudes from our 
spectroscopic analysis, both intrinsic and reddened (for the ob- 
served magnitudes we used those from Udalski et al. (2001) for 
the reasons cited above, but we used ou r own value for the B 
band, as the values from lSandagel (Il97lh are old measurements 
that had to be calibrated through a complicated process (see Sect. 
III.c in Sandage's paper)). We see that the model reproduces the 
observed data well after including extinction, although an ex- 
tra absorption would be needed in all bands relative to the K- 
band. This can be due to a number of factors: a slightly different 
reddening law at low Z, extra emission in the K band, or extra 
absorption at shorter wavelengths, produced by a disk, or a com- 
bination of all of them. 

In Fig. [To] we show the comparison of selected lines in the 
observed spectrum with the CMFGEN model, after being con- 
volved with the appropriate instrumental resolution. Although 
abundances have not been determined, implying that the fit can- 
not be perfect, we see that the general behavior is well repro- 
duced by the model. No effort to determine the rotational veloc- 
ity has been carried out, as the resolution is modest. Parameters 
in Table |2] have been adopted as the final stellar parameters. 

We see that the fit to the observed spectrum is quite good, 
although there is a trend for the absorption part of the higher 
Balmer lines to be broader than the model calculation. This can- 
not be attributed to terminal velocity uncertainties as the lower 
members fit well (note that Galactic stars would have spectra 
redshifted with respect to IC 1613 stars). In a star+disk scenario, 
this could be due to the differential absorption produced by the 
disk in the series members combined with the different levels 
of disk continuum. We also emphasize that the absorption part 
of some Fell P-Cygni profiles (Feii AA4924, 5018) is predicted 
too strong as compared to observations. This is not solved by 
changing the stellar parameters, and thus we rely on the emission 
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peaks. These points might indicate that the density structure of 
V39 is not as simple as described by our spherically symmetric, 
one-dimensional models. The presence of a disk is not the only 
possibility: wind cl umping, for example, was not considered. As 
is well known (see iPuls et akl (12008) ). the presence of clump- 
ing may modify the mass-loss rates derived from the H^^ line, for 
which we obtain a value that is higher than expected for a star 
of the given luminosity and metallicity. Whether this represents 
a challenge for the theory of radiatively driven winds cannot be 
said at the present stage. 

We derive a rather typical terminal velocity for an LB Vc, 390 
km s and thus find no compelling evidence that Vm increases 
with d ecreasing metallicity, as suggested by I zotov & Thuanl 
( |2009|) . Increasing the terminal velocity results in too broad ab- 
sorption in and H g and does not help further. We note that 
iDrissen et al. 1 d2001h also found a modest terminal wind velocity 
for the low-metallicity LBV NGC 2363-Vl (290-325 km s"'). 

We have not determined metal abundances, but we have 
constrained the Fe abundance of V39, which is an impor- 
tant parameter, as Fe is the main driver of the stellar wind. 
We obtained (Fe/Feo)= 0.2, which agrees with the value by 
iTautvaisiene et al."* (2007), who obtained a mean iron abundance 
of [Fe/H]= -0.67+0.09 from the analysis of three M-supergiants 
in IC 1613, which corresponds to 0.21 F cp and the present value 
of Z= 0.004 ([Fe/H]= -0.7) derived bv Skillman e t al. I (l2003h 
from isochrone fitting to the young stellar population of IC 1613. 
All three metallicities are higher than the metallicity derived 
from nebular and stellar O, for which we find in t h e hterature val- 
ues between 7.60 and 7.90 (see lLee et aTl (l2003h . lBresohn et al. I 
(120071) . among other results), which corresponds to 0.08 - 0.16 
times the solar ox ygen abundance (assuming the solar O abun- 
dance from Aspl und et al. I ( l2009h ). The fact that the B super- 
giants and H ii regions give consistent O abundances and that the 
same consistency is found for the Fe abundances from (one) hot 
, (three) cool supergiants and the young stellar population is a 
remarkable result, which supports the conclusion that the [a/Fe] 
ratio in IC 1613 is slightly lower than solar. 

Given the number of parameters and the complexity of the 
models, a proper error study is not possible. Within the parame- 
ter domain the star is placed and for the assumed H/He ratio, the 
errors in the stellar temperature are narrowed down to ±1000K 
given the presence of weak He i 5876 and absence of other He i 
lines. The errors in radius and luminosity are dominated by the 
photometry and reddening uncertainties, while the error in radius 
dominates the error in mass-loss (transformed radius). Because 
of the large number of Fe ii lines in the spectra of V39 the Fe 
abundance can be well constrained once the Tgff and M have 
been determined. Thus, 0. 15 dex can be regarded as a conserva- 
tive error on the Fe abundance. 



U ' B 'V 'R' 1 ' ' ' J ' 'H' ' K : 




LOG (WAVELENGTH) 

Fig. 9. Comparison of observed magnitudes of V39 listed in 
Table [U (triangles) with those of our model, both intrinsic (open 
squares) and reddened (circles). Wavelength is given in A and 
the photometric bands are indicated on top of the plot. See text 
for details. 
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Fig. 10. Fit to the observed spectrum of V39 calculated with 
CMFGEN. Stellar parameters are given in Table |2l 



6. Discussion 

In Fig.[TT]we have plott ed evolutionary tracks including rotation 
for the SMC metallicity dMevnet & Maeder .2005) together with 
the position of V39 from the model parameters and some known 
LBVs. From its position in the diagram we can estimate the ini- 
tial stellar mass to be ~25 Mq. According to the evolutionary 
tracks, the star should be in a redward excursion, without possi- 
ble return to the blue. Because of the low metallicity and conse- 
quently weak wind, the mass-loss in the model up to its present 
stage is very small, and therefore the present mass predicted by 
the models is not very different from the initial mass. 

However, the star may have suffered previous episodes of 
mass loss, and the possible presence of a disk may be an in- 



dication of it. At the derived stellar temperature, the Nai lines 
are expected to be much weaker than observed. The spectral 
type should be at least F 5, according to the stellar atlas of 
iDansk & Dennefeld Id 19941) (that contains Galactic stars) to have 
comparable stellar lines. Therefore we expect the Nai D lines to 
be formed in the circumstellar material, and we conclude that 
their presence is consistent with (but not a proof of!) previous 
episodes of mass loss. 

The derived luminosity is at the lower limit of the known 
LBV and LBV candidates. The absolute visual magnitude we 
derive is Mv= -8.06, which is modest for an LBVc, but again 
within the expected range. The position of V39, close to the 
low luminosity end of known LBVs, adds interest to this ob- 
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Star 


Trmeff 


log (L/Lo) 


R/Ro 


M. 


He/H 


M 


v„ 


(Fe)/(Fe)o 


Reference 




(K) 






(Mo) 




(Mo yr-') 


(km s-') 






V39 


9260 


5.34 


183 


25 


0.4 


3.0x10^" 


390 


0.2 


this work 


HD 160529 


8000 


5.60 


330 


13 




1-3x10^^ 


120 




Sterkenetal. (1991) 


R 110 


7600 


5.46 


310 


10 




1x10^^ 


80 




Stahletal. (1990) 


R40 


8700 


5.61 


280 


16 




8x10"* 


125 




Szeifert et al. (1993) 



Table 3. Synthetic photometry for the model given in Table |2l Values are unreddened. 



Mv 


U 


B 


V 


R I 


J 


H 


K 


-8.06 


15.54 


16.32 


16.21 


16.01 15.91 


15.83 


15.67 


15.47 



ject, as the number of known low-lu minosity LBV s (or L BV 
candidates) is very small. In Fig. 4 of ISmith et al. I (l2004l) we 
see only seven objects in this region of the HRD, including 
both LBVs and LBV candidates, belonging to the Milky Way, 
the LMC and the SMC. In Fig. [TT] we have also marked the 
position occupied by NGC 2363-Vl (Drissen et al. 2001) and 
R127 duri ng outburst, an d the lo w-luminosity LBVs HD 160529 



(Galac tic. ISterken et al. I (Il99lh ). RllO (or HD 269662 in the 



LMC. IStahl et al. I (119901) ) and R40 (or HD 6884 in the SMC, 



ISzeifert et al. I ( 1993h ). While the first two are much more lumi 



nous than V39, the latter ones are closer to it, and we are tempted 
to conclude that they are in a similar stage, although a compari- 
son of the spectrum of V39 with those availa ble in the literature 
for HD 1 60529 dC hentsov et al. 2003), RllO dStahl etal.ll990l) 
or R40 (Szeifert et al. 1993) results in clear differences: the P- 
Cygni profiles of Balmer and Fe lines are more pronounced in 
V39, which could indicate that the star is closer to its Eddington 
limit. This would be consistent with V39 having ejected a large 
amount of matter, indicated by its intense reddening (the redden- 
ing towards HD 160295 is more intense, but this is a Galactic star 
at low Galactic latitude, I- -1.73) and having the largest mass- 
loss rate, in spite of its lower metallicity. Whether this is related 
to the main difference, namely the strong, short-period photo- 
metric variations present in V39, cannot be clarified with the 
present data. However, it is clear that the derived mass-loss rate 
is very high for a star of the given luminosity and metallicity. 
Table |2] gives the main parameters of the stars. Note that for HD 
160529, RllO and R40 values were not obtained through a de- 
tailed spectral fit but with a combination of plane-parallel and 
wind models (details can be found in the cited references). 

Alternatively, V39 could represe nt a high-lumin osity version 
of supergiant B [e] stars (sgB[e], see lLamers et al. I (!l998)). par- 
ticularly if we consider that our luminosity might be slightly 
overestimated (we have considered all K-band emission as stel- 
lar). An argumen t suppo rting this possibility is the indication by 
iMiroshnichenko I (12007 ) that these stars are more luminous at 
lower metallicities (but given the uncertainties in the luminosity 
determination and the low sgB[e] statistics we have to be cau- 
tious with this particular argument). When attributing the cooler 
spectral features to the disk and the main photometric variations 
to disk precession, observed features would be compatible both 
with a LBVc and a sgB[e] nature. The main difference would be 
the lower luminosity expected for sgB [e] stars and its larger in- 
frared excess at long wavelengths (because of dust in the circum- 
stellar material). There are no Spitzer or IRAS data providing 
evidence of large IR excesses for V39, and therefore we prefer 
its classification as an LBV candidate, but without completely 
discarding a possible sgB[e] nature which, on the other hand, 
may simply be different evolutionary stages of the same star. 
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Fig. 11. The position of V39 in the HR Diagram (filled circle). 
For comparison we have also plotted (filled triangles) the posi- 
tions of NGC 2363-Vl and R127 during outburst, as well as the 
low luminosity LBVs HD 160529 (during two stages connected 
by a dotted line), and RllO. Also plotted is the position of the 
SMC LBV R40 (with an open square, to avoid confusion). We 
also marked the region occupied by LBVs during outburst (ver- 
tical box) and the one occupied by LBVs during visual mini- 
mum (slanted strip, both from S mith et al. I (l2004l) ). The tracks 
coiTespond to evolutionary models for the SMC including ro- 
tation (with initial velocities of 300 km s '.'Mevnet & Maederl 
(2005)). Initial track masses are given by the numbers left from 
the ZAMS. 



However, we need some force causing disk precession. This can 
be due either to a non-stable disk or to perturbation from some 
companion. In the first case, the disk should have formed re- 
cently (otherwise friction would stop precession) or even being 
still in formation (the high derived mass-loss rate supports this 
possibility). The second case is not discarded by our discussion 
in Sect. |4] as the companion should only be responsible for the 
small perturbation causing precession, and not for the photomet- 
ric changes. 



7. Conclusions 

We have presented intermediate resolution (R=2050-2150) spec- 
troscopic observations of the variable star V39 in IC 1613 se- 
cured with VIMOS at the VLT. The final spectrum is the combi- 
nation of 19 blue spectra and 10 red spectra and is characterized 
by i) strong P-Cygni profiles in the Balmer lines and Fe n lines; 
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ii) an hybrid character, more consistent with an early type super- 
giant (B-A) in the blue and a late type (G) in the red. All spectral 
features are consistent in radial velocity with the star belonging 
tolC 1613. 

The individual spectra do not show significant radial velocity 
variations in spite of the time span covered, which is longer than 
one photometric period (during which photometric changes are 
larger than one magnitude in B and decrease with wavelength), 
and show only moderate spectral line variations. Therefore we 
have discarded scenarios including a binary system, specially 
those demanding large radial velocit y variations. We also d is- 
carded the scenario proposed by Mantegazza et al. I (l2002h (a 
Galactic W Virginis star plus an IC 1613 red supergiant) be- 
cause it is not consistent with the system belonging to IC1613 
and the observations do not show the spectral characteristics and 
changes of a W Virginis star We propose a scenario consisting 
of an early (B-A) type star surrounded by a hot (~ 5000 K) thick 
disk precessing around it. Although this scenario can potentially 
explain the photometric changes and its wavelength dependence 
as well as the lack of significant spectroscopic changes, details 
have to be worked out (for example, the cause of disk preces- 
sion) and independent evidence (like for example the presence 
of double peaks in spectral lines) has still to be obtained. 

The analysis of the spectrum under the assumption that this 
is a single star indicates a low-metallicity (Fe ~ 0.2 Fe©, from 
the iron P-Cygni lines), hot evolved star with a strong wind. 
The agreement of this result with other spectroscopic Fe o r 
isochrone fitting global metallicities (iTautvaisiene et alH (l2007l) . 
[Skillman et al. (2003)) on the one hand and the stella r and nebu- 
lar O abundances on the other (Breso lin et al. I (l2007h . lLee et al. I 
(12003 )) strongly points to a low [a/Fe] ratio in IC 1613 as com- 
pared to solar values, although more work will be needed to con- 
firm this possibility. 

From its position in the HRD the star lies in the low- 
luminosity edge of LBV and LBV candidates and its progenitor 
may have had a mass of about 25 (which would be close 
to the present mass). The derived mass-loss rate is higher than 
expected for a star of the given luminosity and metallicity, and 
it may represent a challenge to theory. The terminal velocity is 
consistent with other LBVc , and we could not confirm the sug- 
gestion bv llzotov & Thuan] (|2009<) tha t it increases at l ow me tal- 
licities. Our result agrees with that of iDrissen et al. 1 (1200 Ih for 
NGC 2363-Vl. If confirmed, V39 would be the resolved LBV 
candidate located in the most metal poor environment known to 
date. Alternatively, the star could be a high-luminosity version of 
a sgB[e]. This second alternative would imply some significant 
dust formation around V39 (which we could neither confirm nor 
reject) and is not inconsistent with the star being an LBV at some 
stage in its evolution. 

The combination of the low abundances with the high mass- 
loss rate raises the issue of how the winds of evolved stars be- 
have at low metallicities (like those of LBV and sgB[e] stars) and 
whether they are driven by line radiation, pulsations or contin- 
uum radiation playing a stronger role than in Galactic stars. Our 
data are not sufficient to explore this question in detail. Higher 
resolution spectra, allowing us to dissentangle the contributions 
of the star from that of possible circumstellar material, disk or 
faint companion are needed. 

Whether V39 is an LBVc or a new transient phase of mas- 
sive star evolution, an LBV impostor, is an important question. 
If the former, then it will hold the record as the lowest metal- 
licity LBVc known which has important implications for stellar 
evolution (as described in the introduction). If the latter, then un- 
derstanding what this phase might represent has great intrinsic 



interest, in addition to raising the issue of whether or not there 
may be other LBV impostors which are undiscovered due to a 
lack of adequate data, like a light curve. Unravehng its nature is 
therefore of great importance. 

Acknowledgements. The authors thank I. Negueruela and S. Simon-Diaz 
for valuable and inspiring discussions. We also thank the referee, Dr 
Laurent Drissen, for many positive and valuable comments during the revi- 
sion of the manuscript. This project has been supported by Spanish grants 
AYA2007-67456-C02, AYA2008-06166-C03-01 and AYA2008-06166-C03-02 
and was partially funded by the Spanish MICINN under the Consolider- 
Ingenio 2010 Program grant CSD2006-00070: First Science with the GTC 
(http://www.iac.es/consolider-ingenio-gtcl. We acknowledge the use of UVES 
POP data obtained under ESO DDT Program ID 266.D-5655 



References 

Antonello, E., Mantegazza, L., Fugazza, D., Bossi. M. & Covino, S.. 1999, A&A, 
349, 55 

Asplund. M., Grevesse, N., Sauval, Scott, P, 2009, ARA&A, 47, 481 

Bagnulo, S., Jehin, E., Ledoux, C. et al., 2003, Messenger 114, 10 

Bresolin, R, Urbaneja, M.A., Gieren, W., Pietrzyhski, G. & Kudritzki, R.P, 

2007, ApJ, 671,2028 
Chentsov, E.L., Ermakov, S.V., Klochkova, V.G. et al., 2003, A&A, 397, 1035 
Clark, J.S., Larionov, V.M. & Ai'kharov, A., 2005, A&A, 435, 239 
Dansk, A.C. & Dennefeld, M., 1994, PASP 106, 382 
Dolphin, A.E., Saha, A., Skillman, E.D. et al., 2001, ApJ, 550, 554 
Drissen, L., Crowther PA., Smith, L.J. et al., 2001, ApJ, 546, 484 
Lebre, A. & Gillet, D., 1992, A&A, 255, 221 
Garcia, M., Herrero, A., Vicente, B. et al., 2009a, A&A, 502, 1015 
Gonzalez-Delgado, R.M., Perez, E., Tenorio-Tagle, G. et al., 1994, ApJ437, 239 
Hillier, D.J. & Miller, D.L., 1998, ApJ, 496, 407 
Izotov, Y.L, Thuan, T.X., & Guseva, N.G., 2007, ApJ, 671, 1297 
Izotov, Y.I. & Thuan, T.X., 2009, ApJ, 690, 1797 

Karachentsev, I.D., Dolphin, A.E., Geisler, D. et al., 2002, A&A, 383, 125 
Kotak, R. & Vink, J.S., 2006, A&A, 460L, 5 

Lamers, H. J.G.L.M., Zickgraf, F.J., de Winter, D., Houziaux, L. & Zorec, J., 

1998, A&A 340, 117 
Lee, M.G., Freedman, W.L. & Madore, B.F, 1993, ApJ, 417, 553 
Lee, H., Grebel, E.K. & Hodge, PW., 2003, A&A, 401, 141 
Lu, N.Y, Hoffman, G.L., Groff, T, Roos, T. & Lamphier, C, 1993, ApJS, 88, 

383 

Luridiana, V., Peimbert, M. , & Leitherer, C, 1999, ApJ, 527, 110 

Makarova, L.N. & Karachentsev, I.D., 1998, A&AS, 133, 181 

Mantegazza, L., Antonello, E., Fugazza, D., Covino, S. & Israel, G.L., 2002, 

A&A, 388, 861 
Meynet, G. & Maeder, A., 2005, A&A, 429, 581 
Miroshnichenko, A.S., 2007, ApJ 667, 497 
Najarro, R, Hillier, D. J., & Stahl, O. 1997, A&A, 326, 1117 
Najan-o, R, Riger, D.R, Hillier, D.J., Geballe, TR. & Kudritzki, R.P, 2009, ApJ, 

691, 1816 

Onifer, A.J. & Guzik, J.A., in lAUS 250, eds. J. Puis, P Crowther and R. Bresolin, 
p. 83 

Owocki, S.P, Gayley, K.G. & Shaviv, N.J., 2004, ApJ, 616, 525 
Petit, v., Drissen, L. & Crowther PA., 2006, AJ, 132, 1756 
Pietrzyhski, G., Gieren, W., Sozyhski, I. et al., 2006, ApJ, 642, 216 
Puis, J., Vink, J.S., & Najarro, R, 2008, A&ARv, 16, 209 

Pustilnik, S.A., Tepliakova, A.L., Kniazev, A.Y. & Burenkov, A.N., 2008, 

MNRAS, 388, L24 
Sandage, A., 1971, ApJ, 166, 475 

Skillman, E.D., Tolstoy, E., Cole, A.A. et al., 2003, ApJ, 596, 253 
Smith, N., Vink, J.S. & de Koter, A., 2004, ApJ, 615, 475 
Smith, N. & Owocki, S. P, 2006, ApJ, 645, L45 

Stahl, O., Wolf B., Klare, G., Juttner A., Cassatella, A., 1990, A&A, 228, 379 
Sterken, C, Gosset, E., Juttner A., et al., A&A, 247, 383 
Szeifert, Th., Stahl, O., Wolf B. et al., 1993, A&A, 280, 508 
Tautvaisiene, G., Geisler D., Wallerstein, G. et al., 2007, AJ, 134, 2318 
Tolstoy, E., Saha, A., Hoessel, J.G., McQuade, K., 1995, AJ, 1 10, 1640 
Udalski, A., Wyrzykowski, L., Pietrzyhski, G, 2001, AcA 51, 221 
Vink, J.S., 2009, astro-ph 0905.3338 



